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0\ ■ ABSTRACT 

Thermal fluctuations in the magnetization of interstellar grains will produce 
magnetic dipole emission at z/ ^ 100 GHz. We show how to calculate absorption and 
emission from small particles composed of material with magnetic as well as dielectric 
. properties. The Kramers-Kronig relations for a dusty medium are generalized to 

^\ I include the possibility of magnetic grains. 

The magnetic permeability as a function of frequency is discussed for several 



o 
o 



' candidate grain materials. Iron grains, or grains containing iron inclusions, are likely to 

' have the magnetic analogue of a Frohlich resonance in the vicinity of ~ 50 — 100 GHz, 

G^ ■ resulting in a large magnetic dipole absorption cross section. 

We calculate the emission spectra for various interstellar grain candidates. While 
"ordinary" paramagnetic grains or even magnetite grains cannot account for the 
2 ! observed "anomalous" emission from dust in the 14 — 90 GHz range, stronger magnetic 

I dipole emission will result if a fraction of the grain material is ferromagnetic, as could 

be the case given the high Fe content of interstellar dust. The observed emission 
from dust near 90 GHz implies that not more than ~ 5% of interstellar Fe is in the 
^ ' form of metallic iron grains or inclusions (e.g., in "GEMS"). However, we show that 

■ if most interstellar Fe is in a moderately ferromagnetic material, with the magnetic 

properties properly adjusted, it could contribute a substantial fraction of the observed 
14 — 90 GHz emission, perhaps comparable to the contribution from spinning ultrasmall 
dust grains. The two emission mechanisms can be distinguished by measuring the 
emission from dark clouds. 

If ferromagnetic grains consist of a single magnetic domain, and are aligned, the 
magnetic dipole emission will be linearly polarized, with the polarization depending 
strongly on frequency. 



Subject headings: ISM: Atomic Processes, Dust, Radiation; Cosmic Microwave 
Background 
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1. Introduction 

Experiments to map the cosmic background radiation (CBR) have stimulated renewed 
interest in diffuse Galactic emission. Sensitive observations of variations in the microwave sky 
brightness have revealed 14-90 GHz microwave emission which is correlated with 100 //m thermal 
emission from interstellar dust (Kogut et al. 1996; de Oliveira-Costa et al. 1997; Lcitch ct al. 
1997). The origin of this "anomalous" emission has been of great interest. While the observed 
frequency-dependence appears consistent with free-free emission (Kogut et al. 1996; Leitch et al. 
1997), the observed intensities cannot be due to free- free emission from interstellar gas: for gas 
temperature T 10^ K the required spatial variations in emission measure would produce ~ 60 
times more Ha than allowed by observations (Gaustad et al. 1996; Leitch et al. 1997), while 
T ?a 10^ K gas would produce powerful X-ray emission which can be ruled out on energetic grounds 
(Draine & Lazarian 1998). 

It has been proposed that the observed 14-90 GHz emission is rotational emission from the 
population of ultrasmall dust grains believed to be present in the interstellar medium (Draine & 
Lazarian 1998a, b). Such electric dipole emission is a natural consequence of current models of 
interstellar dust, and seems likely to be the mechanism responsible for the observed emission. It is 
nevertheless of interest to investigate whether other emission mechanisms exist which would also 
be capable of producing strong emission at these frequencies. 

Here we discuss an entirely different emission mechanism: thermal fluctuations of the 
magnetization within individual interstellar grains. Such thermal fluctuations will result in 
magnetic dipole radiation. The magnitude of the magnetic fluctuations depends upon the 

magnetic properties of the grain material, and we consider a number of possibilities, including 
"ordinary" paramagnetism as well as superparamagnetism, ferromagnetism, and ferrimagnetism. 

If about 30% of the grain mass is carbonaceous, then Fe and Ni contribute nearly 30% of the 
mass of the remaining grain material. It would not be surprising, then, if some fraction of the 
interstellar grain population were strongly magnetic; metallic iron/nickel, magnetite (Fe304), and 
maghemite (7Fe2 03) are plausible candidates to consider. 

At infrared and optical frequencies (u ^ IQ-*^^ Hz) it has been customary to neglect 
the magnetic properties of the grain material when computing absorption or emission of 
electromagnetic radiation by dust grains; this is an excellent approximation since most materials 
have negligible magnetic response to oscillating magnetic fields at frequencies ^ 10^^ Hz. This 
is because "magnetism" is due to ordering of electron spins, and the maximum frequency for 
electron spins to reorient is the precession frequency of an electron in the local magnetic field (due 
mainly to other electron spins) in the material; these gyrofrequencies do not exceed ~ 20 GHz. In 
contrast, the electron charge distribution in grains can respond strongly to applied electric fields 
at frequencies as large as 10^^ Hz. 

At microwave frequencies materials can respond to both electric and magnetic fields, and both 
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the dielectric constant €{uj) and the magnetic permeability fi{uj) can play a role in the absorption 
and emission of electromagnetic radiation. 

The present study is primarily directed at calculation of thermal emission from magnetic 
dust grains. The complex permeability //(w) obtained here is also relevant to the process of 
grain alignment via magnetic dissipation (Davis & Greenstein 1951), particularly for ultrasmall 
grains undergoing very rapid rotation (Draine & Lazarian 1998a, b; Lazarian & Draine 1998) and 
suprathermally rotating graints (Purcell 1979). 

The optical properties of small particles, including both dielectric and magnetic effects, are 
reviewd in In §^ we generalize Purcell's (1969) application of the Kramers-Kronig relations to 
interstellar grains to include the case of magnetic grains. 

The expected form of the frequency-dependence of the complex permeability fi{uj) is 
discussed in §|5[ In §^ we review the properties of ordinary paramagnetic materials, and in 
we discuss various classes of iron-rich materials, including: single-domain iron grains (§ |7.1| ); 
bulk (multidomain) iron ( ^7.2| ); materials with single-domain Fe inclusions (§[7.3|), including the 
superparamagnetic limit where the inclusions contain less than 10^ Fe atoms each (§[7.4|); 
ferrimagnetic materials such as magnetite and maghemite (§ |7.5D ; and antiferromagnetic materials 

(§PD- 



In §^ we consider several ways in which iron can be present in the interstellar grain population, 
and for each case we calculate the expected thermal emission contributed by magnetic fluctuations. 
Iron grains, if present, would produce strong 50 — 100 GHz emission; this can be used to place 
an upper limit of ~ 5% on the fraction of the interstellar Fe which can be in the form of pure 
iron grains or inclusions. In ^8.2| we discuss the possibility that the Fe could be in a less strongly 
magnetic form (e.g., impure iron), and show that there are plausible magnetic properties for which 
such grains could dominate the 14-90 GHz emission. The polarization of the magnetic dipole 
emission from ferromagnetic grains is discussed in §0. 



We discuss our results in SIO, and summarize in Sll 



2. Optics of Magnetic Grains 

2.1. Spherical Grains 

For grains which are small compared to the wavelength, the cross section for absorption of 
electromagnetic waves can be written as the sum of electric dipole and magnetic dipole cross 
sections (Draine & Lee 1984): 

For nonmagnetic materials (i.e., materials with magnetic permeability ^ = I), the magnetic dipole 
contribution arises from induced eddy currents in the grain, which give it an oscillating magnetic 
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dipole moment which is out of phase with the appUed magnetic field; the absorption cross section 
for this case has been given previously (Landau & Lifshitz 1960; Draine & Lee 1984). For materials 
with a magnetic as well as electric response, the absorption cross section for small spheres of 
radius a can be shown to be 
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where V = Aira'^/S, e = ei + ie2 is the dielectric constant,]^ and jj. = 1 + Attx = ij-i + i/X2 is the 
magnetic permeability. The dipole approximation (|l]) is validQ provided 
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In eq. (Q) the first term is due to the polarization in response to the applied electric field, and 
the second term is due to the oscillating circular magnetization induced by the time-dependent 
displacement current dT)/dt. Similarly, in eq. (^) the first term is due to the magnetization in 
response to the applied magnetic field, while the second term is the heating due to the "eddy 
currents" induced by the time-dependent magnetic field d^/dt. By analogy, we will refer to the 
second term in eq. (§) as being due to "eddy magnetization". Under some circumstances the 
second-order "eddy" contributions can be important - for example, "eddy current" dissipation is 
important for a ^ 0.2 /im graphite grains at 3 x 10^^ Hz (Draine & Lee 1984). However, for 
a ^ 3 X 10^^ cm interstellar grains at ^ 10^^ Hz, (uja/c)'^ < 4 x 10~^ and these higher-order 
terms are generally negligible. The "eddy" terms will therefore be neglected in the remainder of 
this paper. Thus, 



C(f) ^ e2|M + 2|2 

abs 



(5) 



2.2. Ellipsoidal Grains 

Ellipsoidal grains with semiaxes ai < 02 < 03 are characterized by "geometrical factors" or 
"depolarization factors" Li > L2 > which depend only on the grain shape (see Bohren & 
Huffman 1983). Table |^ contains Lj for selected ellipsoidal shapes. 



^ The electrical conductivity a = uje2/4:n. 

^ Eq. neglects "shielding" of the grain interior from the applied magnetic field by induced eddy currents, and 
(^) neglects shielding from the applied electric field by induced time-dependent magnetization. It can be shown that 
such shielding is negligible provided the dipole validity criterion is satisfied. 
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For the moment, we assume e and fi to be isotropic; we will discuss anisotropic /i in §^ below. 
The absorption cross section for incident radiation with the electric vector parallel to semiaxis e, 
and magnetic vector parallel to semiaxis h, is 

£2 
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|1 + L;,(M-1)|2 

where V = (47r/3)aia2a3, and we have neglected the "eddy current" and "eddy magnetization" 
terms. The thermal emission from the grain will be polarized. Interstellar grains undergoing 
superthermal rotation will tend to have their short axis oi || J, where J is the angular momentum. 
If J is in the plane of the sky, then the degree of linear polarization will be 
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|1 + L2(^-1)|2 |1 + L3(/U-1)|2. 

From eq. (^0) it is apparent that (assuming isotropic materials) polarization arises only to the 
extent that both (1) |e — 1| and/or |/i — 1| differ from zero, and (2) the geometrical factors Lj differ 
from one another. 

The contribution of magnetic dipole emission to the polarization from Fe grains will be 
discussed in §P 



3. Kramers-Kronig Relations 

The Kramers-Kronig relations (see Landau & Lifshitz 1960) can be used to relate 
total grain volume to the wavelength integral of extinction (Purcell 1969), or to relate the 
wavelength-dependence of linear and circular polarization produced by interstellar grains 
(Shapiro 1975; Martin 1975a, b). The original discussion by Purcell (1969) assumed nonmagnetic 
grains; here we generalize the discussion to include magnetic grains. Following Purcell, 
consider an electromagnetic plane wave propagating through the dusty interstellar medium: 
E = Re(^oe*^^'~*'^*), H = Re(-ffoe*^^~*'^*)- If Xe(w) and Xm(w) are the electric and magnetic 
susceptibility of the medium, then the plane wave dispersion relation is 

k^C^ =U^{l + A'KX,){l+^T:Xn^ . (11) 

For the interstellar medium we can assume \xe\ ^ 1, \Xm\ ^ 1- Consider grains with number 
density Ugr, volume V , and extinction cross section CextW, where A = 27rc/u;. Then 

ngrCextW ~ ^Im[Xe(w) + Xm('^)] • (12) 
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Similarly, the birefringence of the interstellar medium is proportional to Re[xe('^) + Xmi^)]- The 
linear response functions Xe and Xm must separately satisfy the Kramers-Kronig relation, 

2 ujduj 

Relxiiioo)] = -V :Am.[xi{uj)] (for i = e,m) . (13) 

vr Jo ^ - 

where V indicates that the principal value is to be taken. Taking = one obtains the 
generalization of Purcell's (1969) result: 

c 

Re [Xe(0) + Xm(0)] = ^j^ dX ngrCextiX) ■ (14) 

For randomly-oriented ellipsoids with semiaxes oi, 02, 03, we have 

W 

Xe{0) = —ngrF{eo;ai/a3,a2/a3) (15) 
47r 

3V 

Xm(0) = —ng,.F(;Uo; 01/03,02/03) (16) 
47r 

where 

and the Lj are the geometric factors discussed above (see Table |T[). For spheres, Lj = 1/3 and 
F{x) ~ 1 for X ^ 1 [Fig. 1 of Purcell (1969) shows F{x) for spheroidal grains (oi = 02 or 
02 = 03)]. 

Thus 

1 r°° 

ngrV [F(eo) + F(^o)] = ^ Ugr CeAX) , (18) 

where eo and /io are the zero-frequency dielectric function and magnetic permeability of the grain 
material. For nonmagnetic materials, F{fj,Q = 1) = and we recover Purcell's (1969) result. 

Grain materials of interest have eo ^ 3 (conducting grains have eo — > 00). For magnetic grains 
with fiQ ^ 3, F{ij,o) « F(eo), and we see from eq. ( |T8| ) that / CextdX must be larger by a factor 
[1 -|- F / F (eo)] than it would have been had the grain been nonmagnetic. 

Purcell (1979) used / UgrCextdX to estimate UgrV. At first sight eq. (^) appears to reduce 
the required volume of interstellar dust by about a factor of two, but it must be recognized that at 
the wavelengths A < 1 mm which contributed to Purcell's estimate for / CextdX, magnetic effects 
are indeed unimportant, so that 

A second consequence, however, is that we see that if -F(^o) ~ -^(^o); then magnetic dipole 
absorption must contribute a cross section C^f' such that 

ct'Ux ^ I &:^^dX ; (19) 



magnetic dipole absorption and emission must dominate the emission from (stationary) interstellar 
dust grains at the frequencies ^ 30 GHz where the grain material has a strong magnetic response. 
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4. Candidate Materials 

In addition to "ordinary" paramagnetism, grains may exhibit strong magnetic ordering: either 
ferromagnetic, ferrimagnetic, or antiferromagnetic (see Morrish 1980, and below). 

Ordinary paramagnetism (§ |6.lD arises when the interaction producing ahgnment of electron 
spins is just the interaction of the electron's magnetic moment with the local magnetic field. 
Ferromagnetism, ferrimagnetism, and antiferromagnetism (jQ) occur when the electron spins are 
ordered due to the "exchange interaction". By mass, Fe is the fifth most abundant element, 
following H, He, O, and C, and nearly all of the interstellar Fe is in grains (Savage & Sembach 
1996). Any strongly magnetic grain materials will almost certainly contain Fe. Table ^ lists 
a number of possible Fe-containing materials which are magnetically ordered at the T ^ 20K 
temperatures of interstellar dust. 



Equations (|6|,0) show that the absorption cross sections are proportional to the grain volume. 
If 100% of solar abundances (Anders & Grevesse 1989) of Mg, Si, and Fe are incorporated into 
silicates with the approximate composition Mgi.iFeo.9Si04 and a density p ~ 4gcm~^, the 
resulting volume of silicate per H atom is Vq = 2.5 x 10^^'' cm^/H, which we will use as a fiducial 
grain volume per H atom. 

Suppose a fraction /y of the Fe is incorporated into a material Y with density pY, with Fe 
contributing a fraction zy of the mass. Then the volume per H of material Y is 

Vv „ / 4gcm^'^\ , , 

-f = 0.30/y ^ . (20) 

Table | shows Vy jV^ for various candidate minerals. We see that Fe-rich materials are 
likely to make up a significant fraction of the grain volume, and these materials tend to be 
magnetically-ordered. 



5. The Magnetic Susceptibility 

Solids exhibit a variety of different responses to weak applied oscillating magnetic fields. 
When a static magnetic field is present, or the substance is spontaneously magnetized, the 
susceptibility xW) ^ tensor nature (see, e.g., Jones & Spitzer 1967, §IIIa). In the present 
discussion we will treat x ^-s a scalar, but will extend the discussion to spontaneously-magnetized 



material in §7.1 and Appendix 



Aside from the diamagnetism of superconductors and materials with no unpaired electron 
spins, at low frequencies most materials are characterized by x(0) > 0) resulting from either 
"normal paramagnetism" or magnetic ordering in the form of superparamagnetism, ferrimagnetism, 
or ferromagnetism. We consider these cases below, but first discuss the likely behavior of the 
frequency-dependent susceptibility x{^) = Xi + ^X2- 
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ai : 02 


: 03 


Li 


L2 




Note 


1 : 1 


1 


0.3333333 


0.3333333 


0.3333333 


sphere 


1:1: 


L.25 


0.3620042 


0.3620042 


0.2759916 


prolate spheroid 


1:1: 


1.5 


0.3835093 


0.3835093 


0.2329815 


prolate spheroid 


1 : 1 


2 


0.4132180 


0.4132180 


0.1735640 


prolate spheroid 


1 : 1.25 


1.25 


0.3944403 


0.3027798 


0.3027798 


oblate spheroid 


1 : 1.25 


: 1.5 


0.4189528 


0.3233250 


0.2577222 




1 : 1.5 


1.5 


0.4459056 


0.2770472 


0.2770472 


oblate spheroid 


1 : 1.5 


: 2 


0.4837282 


0.3050063 


0.2112656 




1 : 2 


2 


0.5272003 


0.2363999 


0.2363999 


oblate spheroid 



Table 1: Geometrical factors Lj for ellipsoids with semiaxes Oj. 



Mineral Y 


Form*^ 


47rM,(G)'= 


/)y(gcm ^) 




Fe2Si04 (fayalite) 


antiferromag.'' 





4.39 


0.51/y 


FeSiOs (pyroxene) 


antiferromag.*^ 









FeS (troilite) 


antiferromag.*^ 





4.83 


0.39/y 


FeCOs (siderite) 


antiferromag.*^ 





3.96 


0.63/y 


FeMgSi04 (olivine) 


antiferromag.*^ 





4. 


0.92/y 


aFe203 (hematite) 


antiferromag.'' 





5.26 


0.33/y 


FeO (wustite) 


antiferromag? 





5.75 


0.27/y 


Fei_i.S (pyrrhotite) 


ferrimagnetic 


1130^ 


4.6 


0.41/y 


MgFe204 (magnesioferrite) 


ferrimagnetic 


1760 


4.18 


0.51/y 


NiFe2 04 (trevorite) 


ferrimagnetic 


3770 


5.35 


0.47/yg 


7Fe2 03 (maghemite) 


ferrimagnetic 


4780'" 


4.88 


0.36/y 


Fe304 (magnetite) 


ferrimagnetic 


6400 


5.2 


0.32/y 


Fe (iron) 


ferromagnetic 


22000 


7.86 


0.15/y 



Table 2: Magnetic compounds containing Fe. 



"Data selected from review by Carmichael (1989). 
''Magnetic ordering at T < 20K for ideal mineral. 
'^Ms is the spontaneous magnetization at T < 20 K. 

'^Volume W per H atom of grain material Y containing a fraction fy of the solar abundance of Fe, relative to 
Vo = 2.5 X 10"^^ cm^/H, the volume of silicate which can be made using all the Si, Mg, and Fe from solar abundances. 
■^Small particles may be ferrimagnetic (see §7.(3). 
■'^For X = 1/8, where Ms is maximized. 
^ fv < 0.11 required by cosmic abundances. 
''Dunlog & Ozdemir 1997). 
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5.1. Drude Susceptibility 

Suppose that the magnetization M{t) obeys the equation of motion 

M = Lo^o[x{0)H -M]-M/to , (21) 

with three parameters: the static response x(0)) ^ resonance frequency ujq, and a characteristic 
damping time tq. Then we obtain the Drude form for the susceptibihtyj^ 

XM= . , . , (22) 

where we have defined 

T^ii^lror^ . (23) 
The real and imaginary parts of x('-^) s-^'s 

X1=X(0) , /7ty, ,2 (24) 

X2 = X(0) ,- , I ^Ji2 , / X2 • (25) 



A susceptibihty of the form (22) arises, for example, as a solution to the Bloch equations (c.f. 
Pake 1962, eq. 6-24) for a magnetized material. 

The product {jJqt = I/luoTq determines the shape of x('^)- For ll/qt < 2, x(^) has a resonance 
near ujq. For wqt > 2, the system is "overdamped" , and responds to a step function change in the 
applied H with two distinct relaxation times. 



5.2. "Critically-Damped" Susceptibility 

For ujqt = 2, xi^) is "critically-damped", with a single relaxation time ujq^ = 2tq. We will 
consider this as a plausible form for the response function x{^) ^or a system which is essentially 
nonresonant. For this case x takes the simple form, shown as the heavy curves in Fig. |l|: 

x"''' = 7r^^ ' (26) 

(1 — liOT/Zj^ 

' ^ ^ [1 + (a;T/2)2]2 V ; 



^Morrish (1980) suggests a Lorentzian form for the susceptibility, corresponding to an equation of motion 
M = [x{0)H — M]/t. Because this lacks an "inertial" term, we believe that it overestimates the magnetic response 
at very high frequencies. 
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Fig. 1. — Functional form considered for the magnetic susceptibility X = Xi + ^X2) for different 
values of the parameter uot [see eq. (22 - 25)]. 
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x(o)- 



LOT 



(28) 



[1 + (W2)2]' 

The superscript "(cd)" stands for "critically-damped". We will use x^'^^'^ to estimate the response 
associated with the magnetization of a single-domain particle, or of paramagnetic materials. We 
will identify the characteristic frequency loq = 1/2tq with the gyrofrequency of an electron in the 
internal field. From eq. we see that the response falls off rapidly when the frequency 

to > LUq = 2/t. 



6. Paramagnetism 
6.1. Ordinary Paramagnetism 

For ordinary paramagnetism, the zero-frequency susceptibility is (Draine 1996) 

X(0)«4xl0-Vp(^)'(^) , (29) 

where fp is the fraction of the atoms which are paramagnetic, with magnetic moments pub, 
where the Bohr magneton fiB = eTi/2m(,c. If essentially all of the interstellar Mg, Fe, and Si 
are incorporated into material with approximately the composition of MgFeSi04, then a fraction 
fp ~ 1/7 of the atoms would be Fe, presumably in the form of Fe^'''(^i^4) or Fe^^ (^6*5/2) ions; 
these have p ~ 5.4 and 5.9, respectively (Morrish 1980). Fayalite (Fe2Si04) is antiferromagnetic 
(Carmichael 1989) and it appears that intermediate olivines (Mg2;Fe2-a;Si04) will also be 
antiferromagnetic (Duff 1968). The magnetic character of amorphous olivine is uncertain. While 
weak ferrimagnetism seems likely (see below), here we consider paramagnetic behavior. 

The damping time r is the spin-spin relaxation time, which is essentially the time for electron 
precession in the random magnetic fields within the solid (see Caspers 1967). For amorphous 
olivine we estimate the r.m.s. gyrofrequency to be (Draine 1996) 

u;o/2vrf«8GHz . (30) 

Draine (1996) used a different functional form|^ to estimate the paramagnetic susceptibiity, but 
we now consider the "critically-damped" form ;^('^'^) [eq.(|26|)] to provide a better estimate at 
frequencies uj » ujq. In Fig. |2| we show ^2 = ^t^X2'^\ for = 2/loq = 4 x 10^^-*^ s. 

At frequencies v ^ 200GHz "astronomical silicate" (Draine & Lee 1984) has ei ~ 11.6, 
€2 ~ 1.13 X 10^^(i//GIIz). From eq. (|), it can then be seen that C^Xf' > C^^^f when fi2 > •05e2 
(assuming fii ~ 1, as is the case for ordinary paramagnetism), so for reference we also show .05e2 
m Fig. |. From Fig. | it can then be seen that if FeMgSi04 grains are merely paramagnetic, the 
thermal emission from these grains at 30 GHz will evidently be dominated by magnetic dipole 
emission arising from thermal fluctuations in the magnetization. 



^ The form used by Draine (1996) had X2 oc a; ^ for uj 00. 
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0.01 0.1 1 10 102 

i/(GHz) 

Fig. 2. — Estimated imaginary component of magnetic permeability for "normal" paramagnetism 
in silicate grains (curve labelled PM). The dot-dash line labelled ".05e2" shows the level above 
which magnetic dipole absorption due to paramagnetism dominates electric dipole absorption (see 
eq-i- 
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6.2. Stark Splitting in Paramagnetic Materials 

Additional magnetic absorption can arise in paramagnetic material when crystalline electric 
fields produce Stark-effect splitting of the magnetic sublevels, and this splitting corresponds to 
microwave frequencies (Kittel &: Luttinger 1948). Electric dipole transitions are forbidden, since 
the sublevels have the same parity, but magnetic dipole transitions are allowed. 

This type of absorption has been studied in the so-called Zero Field Resonance experiments 
(see Bramley & Strach 1983 and references therein). For Fe^^ resonances were observed at 24.2 
and 35.4 GHz in natural amethyst and at 7.1 and 8.8 GHz for synthetic brown quartz. Similar 
measurements in sapphire revealed resonances at 11.8 and 18.9 GHz. 

For pure single crystals the absorption is localized in lines, while a continuum of absorption 
features is expected for impure and noncrystalline interstellar grains. The strength of this 
absorption can be estimated (see Appendix At the temperatures of interstellar grains, it 
appears that a permeability fi2 ^ 10~^ would result. We see from Fig. ^ that this would be a 
minor contribution to grain absorption and emissions compared even with what we estimate for 
"ordinary" paramagnetism at 1 — 30 GHz, and far below what we estimate in §0 for ferromagnetic 
or ferrimagnetic materials at microwave frequencies. Accordingly, as a continuum emission process 
this mechanism appears to be unimportant. 

The absorption due to Stark splitting could be concentrated in a few absorption features, 
rather than being spread into a continuum. Even so, the effect will be weak. If, say, 5% of the 
atoms had Stark splitting features within a frequency range Alnw = 0.01, then we see from eq. 
(A4) that fi2 ~ 0.14(|/i„ip//i^) within the narrow feature. At, say, 5 GHz this is only a few times 



larger than the value of ij,2 which we estimate for ordinary paramagnetism (see Fig. |2|) and hence 
would result in a narrow emission feature with central intensity only a few times what we estimate 
for the continuum due to ordinary paramagnetism, which we shall see below is itself two orders 
of magnitude weaker than the observed continuum emission. Nevertheless we cannot exclude the 
possibility that future sensitive measurements of the microwave background might disclose weak 
spectral features arising from the Stark effect. 



6.3. Magnetic Materials in Interstellar Grains 

As noted in the substantial fraction of Fe in interstellar grains creates the possibility 
that some fraction of the grains could be magnetic, with the "exchange interaction" resulting in 
ordering of the atomic magnetic moments within a single magnetic domain (see Morrish (1980) for 
a review of the physics of magnetism, and Dunlop and Ozdemir (1997) for an excellent review of 
magnetic minerals). 

Ferromagnetic materials as components of interstellar grains were apparently first considered 
by Spitzer &; Tukey (1951), who discussed the possible formation of ferromagnetic materials 
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in interstellar grains as the result of grain-grain collisions: metallic iron, Fe304, 7Fe203, and 
MgFe204 were considered likely products. Jones & Spitzer (1967) suggested that interstellar 
grains might contain very small clusters of magnetic materials such as Fe304 or 7Fe203. Shapiro 
(1975) proposed that the observed polarization of starlight could be produced by platelets of 
Fe304. Sorrell (1994) presented a model for the origin of small Fe304 clusters in H2O ice mantles 
irradiated by cosmic rays. Ferromagnetic properties of mixed MgO-FeO-SiO grains were discussed 
by Duley (1978). 

Bradley (1994) argued that certain interplanetary dust particles ("Glass with Embedded 
Metals and Sulfides", or "GEMS") consisted of aggregates of ~ 0.1 fiui interstellar grains, with 
nanometer-sized Fe-Ni metal inclusions. Martin (1995) argued that the properties of these 
grains were consistent with being interstellar dust, and that the Fe-Ni inclusions could be 
superparamagnetic and able to bring about alignment with the galactic magnetic field. Goodman 
&: Whittet (1995) noted that the numbers of such inclusions were consistent with the requirements 
of Mathis' (1986) hypothesis to account for the dependence of degree of alignment on grain size. 

The S in GEMS appears to be in the form of FeS; if the excess of Fe over S is metallic, then 
the atomic abundances given by Bradley (1994) for 3 bulk GEMS indicates a volume filling factor 
(j) ~ 0.03 for metallic Fe-Ni; we will use this value when discussing the properties of silicate grains 
with metallic inclusions (§ [7.31) . 

7. Ferromagnetic and Ferrimagnetic Materials 

There are three classes of magnetically-ordered materials. In ferromagnetic materials the 
atomic spins within a domain are parallel; in ferrimagnetic and antiferromagnetic materials 
the magnetic ions are located on two magnetic sublattices of oppositely directed spins and 
magnetic moments. In antiferromagnetic materials the magnetic moments of the two sublattices 
are exactly opposite in direction and magnitude, so the net magnetization is zero (for a perfect 
crystal at zero temperature). In ferrimagnetic materials the magnetic sublattices do not exactly 
compensate, either because the magnetic moments of the two sublattices differ in magnitude 
(normal ferrimagnetism) or are not precisely opposite in direction (spin-canted ferrimagnetism) . 

7.1. Single-Domain Iron 

Applied magnetic fields change the magnetization of a bulk sample by two processes: rotation 
of the magnetization within a domain, and motion of domain walls. Grains smaller than a critical 
radius Oc will contain only a single domain; for Fe Oc ~ 3 x 10~^ cm (Morrish 1980). Such 
single-domain behavior can take place for clusters of materials which are ferromagnetic (e.g., Fe) 
or ferrimagnetic (e.g., Fe304, or 7Fe203). Spontaneous magnetization has been observed for Fe 
clusters as small as Nd = 20 atoms (Billas, Chatelain & de Heer 1994) 
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If the magnetic material consists of a single domain, then there is clearly no motion of domain 
walls involved in the magnetization process, unlike the case for bulk samples, so susceptibilities 
measured for bulk samples are inapplicable. 

When a ferromagnetic or ferrimagnetic material consists of a single domain, it will be 
spontaneously magnetized with magnetization Mg along one of the "easy" directions. If a 
weak oscillating field H is now applied, the single-domain sample will exhibit a susceptibility 
which depends on the direction of H. In Appendix ^ we show that at low temperatures only 
the component of H perpendicular to the spontaneous magnetization induces a change in 
magnetization of the sample: AM = x±H_l. For Fe, we estimate x±(0) ~ 0.3xbuik{0) (see 
Appendix 

Since changing the magnetization of the sample involves only small reorientations of 
the magnetic moments, and no motions of domain walls, we assume that the only frequency 
characterizing the response is the gyrofrequency ujg = (e/mec)(47rMs/3) of a magnetic moment 
with 5( ~ 2 in the internal field 4:TtMs/3 to which each of the dipoles is subject. For Fe, with 
inMs = 22kG, we estimate a;o/27r ^ 20 GHz. We use x''"'^^ with r ^ 2/wo ~ 1-6 x 10"^^ s 

to estimate x±{^) (see Fig. |^. 

The energy dissipation rate is proportional to lm{fj,±)H'j_. For isotropic radiation incident on 
the grain, {H'j_) = (2/3) (ff^). Assuming the higher order term in eq. (^) to be negligible, we would 
therefore take the angle-averaged magnetic dipole absorption cross section for randomly-oriented 
spheres to be 

whereas for randomly-oriented ellipsoids with semiaxes ai < 02 < as, magnetized along the long 
axis, 

(crj^^E-TT^ r^^i^^-TT^-r-r (32) 



3c 



with IJ.± = 1 + 4:7rx± obtained using eq. (p^) with either x±(0) obtained from eq. (| 



7.1.1. Frohlich Resonance Condition 

We now observe that if our estimate for x±{^) iii Fig. |^ is correct, then Fe particles will have 
the magnetic analogue of a Frohlich resonance (Frohlich 1949; Bohren & Huffman 1983) - a peak 
in Cabs where fi±^i = 1 — LJ^, thus minimizing the denominator in eq. (|^ or in eq. (|3^). 

It is apparent that for the functional forms in Fig. [l| there are either zero or two frequencies 
where this condition is satisfied; if there are two, then Cabs peaks at the higher of the two 
frequencies, which we denote uip. 

The Frohlich resonance condition for spheres requires that fii < —2 for some range of 
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Frequency co/2n (GHz) 

Fig. 3. — Real and Imaginary parts of the permeability /x = 1 + 4:Trx± for single-domain Fe 
with H perpendicular to the direction of spontaneous magnetization (see text). The estimated 
gyrofrequency LOg (where iJ,±^i = 1) and Prohlich frequency ujf (a*±,i = —2) are indicated. 
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frequencies, or xi < —S/Att^ Whether or not Fe particles will have a Frohlich resonance will 
depend on the details of the frequency-dependence of which at this time is not experimentally- 
determined. We are therefore forced to rely on model-dependent estimates. With our estimate 
of X± (see Fig. we estimate that, for spheres, ujp /27r w 70 GHz. For ellipsoids with moderate 
axial ratios, the resonance will be shifted slightly. For example, for a 1:1.5:2 ellipsoid, magnetized 
along the long axis, the Frohlich resonances along axes 1 and 2 occur for fii = —1.07 and —2.28, 
or ujf/2tt^ 100 GHz and 60 GHz. 

Recalling the Kramers-Kronig relation (Landau & Lifshitz 1960) 

Xi(w) = -V / :rdx (33) 

TT Jo — o;^ 

(where V denotes principal value integral) and the fact that X2 > 0, one sees that the integrand 
is negative for x < u;. Quite generally, X2('^) will have X2(0) = 0, will increase with increasing uj 
(linearly at low frequencies), peak, and then eventually decline to zero at high frequencies. If X2 
declines more rapidly than w"^ at high frequencies, then we expect to have xi < above some 
frequency. Because metallic Fe has xo ^ it appears likely that there will be a frequency range 
where 47rxi < —3, so that the Frohlich resonance condition can be satisfied. 

The functional form x^'^'^^ used here provides what we consider a reasonable estimate for x±) 
but the predicted Frohlich frequency ujf/2tt ~ 70 GHz is model-dependent and should be regarded 
as quite uncertain. 



7.2. Bulk Iron 

The static susceptibility of bulk, metallic Fe is x(0) ~ 12. The susceptibility of Fe as a 
function of frequency has been measured using fine wires, thin lamina, and powders (Allanson 
1945; Epstein 1954) up to ~ 10 GHz. The response of bulk material will result from motion 
of domain walls as well as the deflections of magnetization within single domains as discussed 
above. We treat these two responses as additive: x(bulkFe) = (2/3)x± + Xdw, where Xdw is the 
contribution from motion of domain walls. 

The experimental results show considerable scatter, but we can reproduce the general trends 
if Xdw is given by eq.(^) with x(0) ~ 10, r 1 x 10~^ s~^, and ~ 10^*^ s~^ In Fig. |^ we plot //2 
estimated for multidomain metallic Fe using these parameters. 



5 For x'''"^H^) [eq- (||)] it is easy to show that the Frohlich resonance condition for spheres, fii = —2, can be 
satisfied provided x(0) > B/tt = 1.91, with up > \/Sujo. However, up is associated with a distinct peak in Cabs only 
for x(0) ^ 3. 
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Frequency cj/Stt (GHz) 



Fig. 4. — Real and imaginary permeability /x for multidomain Fe (see text). The gyrofrequency 
LOg and Frohlich frequency uf are shown. 
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7.3. Fe Inclusions 

Consider spherical single-domain ferromagnetic inclusions with volume filling factor 
(p distributed in a nonmagnetic matrix. Suppose each randomly-oriented inclusion to be 
spontaneously magnetized along one of its "easy" directions. We may suppose that, in effect, 
(2/3) of these inclusions have their magnetization M perpendicular to the direction of the applied 
H, with the remaining 1/3 either parallel or antiparallel to H (and therefore not contributing to 
the susceptibility). 

The effective susceptibility of the composite material may be estimated using effective medium 
theory (Bohren & Huffman 1984). For the present case of spherical inclusions with small volume 
filling factor (p <^ 1, Maxwell- Garnett effective medium theory is appropriate, so we estimate 

^ (2/3)0X± ,o4^ 

l + (47r/3)x±(l-2,/./3) ^ ^ 

where x± is given by x^'^'^^ [eq.(]2^)]. In Fig. ^ we show ^ estimated for Fe inclusions with a volume 
filling fraction (f> = 0.03 (the estimated filling fraction of Fe-Ni inclusions in GEMS; see § |6.3| in a 
nonmagnetic (e.g., silicate) medium. Note the peak in the effective fi2 for the composite medium 
at uj/2tt ~ 70 GHz; this arises because of the Frohlich resonance in the individual single-domain 
Fe inclusions. 



7.4. Superparamagnetism 

If the Fe atoms are aggregated into single-domain clusters which are sufficiently small, 
then thermal fluctuations will cause the magnetization to fluctuate in direction, with the most 
probable magnetization directions being the most energetically favorable ones. If the inclusions 
are sufficiently small {N^i ^ 6 x 10^ for Fe at room temperature), the energy barrier associated 
with reorientation of the magnetization is small enough that thermal fluctuations can reorient 
the magnetization of a single particle on time scales of seconds or less. Such clusters are termed 
superparamagnetic. If a field H is applied, the change in magnetization will be greater than would 
be the case for larger clusters, where only a small deflection in the direction of magnetization takes 



place, as discussed in §7.1 



Suppose that the grain consists of a nonmagnetic (or weakly magnetic) material containing 
superparamagnetic clusters of Nd Fe atoms, with volume filling factor (psp- A single grain then 
contains N ^ 3.5 x 10^(j)spN^i^{a/10'^ cm)^ clusters. 

The susceptibility of the medium may be approximated as 

^^l + (4vr/3)x±(O)[l-(2/3)0,p]^^^^ ' ^ ^ 

The first term is what would be expected if the individual clusters have their magnetization 
locked into the "easy" direction, with a small reorientation of the magnetization resulting from 
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Fig. 5. — Same as Fig. ^ but for material with 3% of the volume due to inclusions consisting of 
Nci Fe atoms. The Frohlich frequency cop is that estimated for single-domain Fe (see text). 
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the transverse component of the apphed H [see eq. (^)]. The second term is what is expected 
from ahgnment of individual clusters assuming them to be free to orient their magnetization in 
any direction. We note that some previous discussions of the susceptibility of superparamagnetic 
materials (e.g., Jones &: Spitzer 1967; Lazarian 1995; Draine 1996) have omitted the first term 
when considering superparamagnetism. In principle, these two estimates should not be treated 
as simply additive, but in general one or the other will dominate, and factor-of-two accuracy is 
sufficient at this time. 

The superparamagnetic contribution at zero frequency is estimated to be 

Xsp{0) ~ .O350,piV,/ (^^) . (36) 

In general (even if spherical), a cluster must overcome an energy barrier in order to substantially 
reorient its magnetization. Laboratory experiments show the relaxation process to be be thermally 
activated, with a characteristic relaxation rate 

^ Aexp{-N,i9/T) , (37) 

with A ~ 10^ s~^ and 9 ~ O.OllX for metallic Fe spheres (Bean & Livingston 1959; Jacobs & 
Bean 1963). We will employ ( p2|) with ujq = lO^'^s"^ to estimate Xsp{^)- 

Fig. P shows fi2 for silicate grains containing superparamagnetic Fe clusters with volume 
filling factor (psp = 0.03. Such inclusions, if present in all of the silicate grains, would account 
for / = (/)sp/0.15 = 20% of the total Fe abundance (see Table ^). Permeabilities are shown for 
Nci = 10^, 10^, and 10^ Fe atoms per cluster. 

The superparamagnetic response is large at low frequencies, but is small compared to normal 
paramagnetism at frequencies ^ 1 GHz, because the estimated relaxation rate is too slow; at 
frequencies ^ 1 GHz, the magnetic response is essentially that of the individual domains having 
their spontaneous magnetizations deflected slightly by the transverse component of the magnetic 
field. 

We conclude that while superparamagnetism may lead to enhancements of magnetic 
dissipation at z/ ^ 1 GHz, at higher frequencies it appears unlikely to make an appreciable 
contribution to absorption by interstellar grains. Therefore the high-frequency response of grains 
with Fe inclusions depends only on the volume fraction (pgp of Fe inclusions, but is insensitive to 
the number N^i of Fe atoms per inclusion, provided only that N^i ^ 20 so that the inclusions are 
ferromagnetic. 



7.5. Ferrimagnetic Minerals: Magnetite and Maghemite 

Magnetite (Fe304) is a commonly-occuring terrestrial mineral, and is a plausible interstellar 
grain material. The low temperature spontaneous magnetization Mg ~ 480 G, from which we 
estimate a gyrofrequency ujg /2tt 6 GHz. 
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At room temperature magnetite is cubic, and the "easy" direction is (111). At ~ 120 K it 
undergoes a phase transition with a shght distortion of the unit cell from cubic to monoclinic 
symmetry; near this temperature the magnetic anisotropy parameter Ki (see Appendix |B|) changes 
sign and the "easy" direction changes to (100). The magnetic anisotropy parameters have been 
measured down to 77 K (Syono & Ishikawa 1963). Extrapolation to T ~ 18K is uncertain; we take 
Ki ~ 1.5 X 10^ ergscm^^, from which we estimate x±(0) ~ 0-83 using eq. (|B4|). 

We have not found data on the high frequency susceptibility of magnetite. Taking 
T = 2/ujg = 5 X 10^^^ s, we estimate /i for single-domain Fe304 using ^(^'^) [eq. (p6|)]; the results 
are shown in Fig. |6|. 

The low temperature oxidation product of magnetite is maghemite, 7Fe203, for which we 
estimate ujg/2ir ~ 4 GHz. From Mg and the crystalline magnetic anisotropy Ki we would estimate 
X±(0) ~ 0-4 (see Appendix |B|). However, measurements by Valstyn, Hanton, & Morrish (1962) 
of the permeability as a function of frequency are better reproduced by a slightly larger value, 
X±(0) ~ 0.6, which we adopt. Our estimate for ^{uj) for 7Fe2 03 is shown in Fig. ^. Valstyn et al. 
measured fi{uj) for a powder of 0.03 — 0.2 ^m diameter 7Fe2 03 spheres suspended in paraffin wax 
with a volume filling factor 0.15; thus in effect there was a volume filling factor of 0.10 for spheres 
with spontaneous magnetization perpendicular to the applied oscillating field. In Fig. ^ we show 
their experimental results for /x — 1 divided by 0.10. The model which we are using appears to 
be in reasonably good agreement with the measured frequency-dependent susceptibility for this 
material. 



7.6. Antiferromagnetic Materials 

In antiferromagnetic materials the exchange interaction leads to magnetic domains with 
magnetic moments which are ordered, but in such a way that the net magnetization in the domain 
is zero. The olivine fayalite (Fe2Si04) is an example of an antiferromagnetic material. These 
substances have zero frequency susceptibilities x(0) which are similar to (but somewhat smaller 
than) those of normal paramagnetic substances with similar concentrations of magnetic ions. Our 
estimate for "normal paramagnetism" would therefore be a reasonable guide to antiferromagnetic 
substances such as fayalite. 

However, we do not expect interstellar grains to contain perfect crystals; defects and impurities 
seem likely to be common, and the materials may be amorphous (as appears to be the case for 
interstellar silicates, as indicated by the profile of the interstellar 10 /im absorption feature due 
to the Si-0 stretching mode). Furthermore, if the magnetic material is in very small inclusions, 
perfect pairing of spins will not be possible along the boundary of the inclusion. For instance, 
Schuele & Deetscreek (1962) find that small (^ 10~^^ cm) particles of NiO are weakly ferro- 
or /ernmagnetic although bulk NiO is anii/erromagnetic. Antiferromagnetism therefore seems 
unlikely - we instead expect amorphous fayalite, for example, to be at least weakly ferrimagnetic. 
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Fig. 6. — Same as Fig. |3| but for single-domain magnetite (Fe304), an example of a ferrimagnetic 
material. Also shown are experimental results for 7Fe203 from Valsytn et al. (1962). 
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The high frequency response of such small inclusions would be determined by the single-domain 
susceptibility ( ^7.3| ). 



For a nonrotating grain whose internal degrees of freedom are in thermal equilibrium at a 
temperature T, the power radiated in frequency interval di' is simply 



where Cabsit^) is the absorption cross section at frequency u, and By{T) is the Planck function. 

We will use the term "vibrational emission" to refer to emission arising from thermal 
fluctuations in the charge distribution, and therefore the electric polarization, of the grain. This 
will produce electric dipole emission. 

We will use the term "magnetic dipole emission" to refer to emission arising from thermal 
fluctuations in the magnetization of the grain material. Estimation of this emission is one of the 
goals of the present paper. 

We use the term "rotational emission" to refer to the emission arising from rotation of the 
grain. This is primarily due to the electric dipole moment which rotates with the spinning grain 
(Draine & Lazarian 1998a, b). In principle, there will be rotational emission arising from the 
rotating magnetic dipole moment of a spinning magnetized grain. Only a ^ 10~^ cm grains spin 
at frequencies v ^ 10 GHz. The magnetic moment of a spontaneously magnetized grain is just 
MsV = (l/3)(47rMs/103G)(a/10-'^cm)3debye. The largest value of Ms is for metallic Fe (see 
Table an Fe grain would have a magnetic moment MgV = 7.3(a/10~^)^ debye, somewhat 
smaller than the electric dipole moment 9.3(a/10~^)^/^ debye estimated for neutral grains by 
Draine & Lazarian (1998b). Other magnetic materials have smaller values of Mg, and hence the 
rotational emission at 10 GHz is expected to be dominated by electric dipole radiation. 



We consider 4 possible components of the interstellar grain population: 
1. 100% of the Si and Fe incorporated into amorphous silicate grains with paramagnetic 



8. Microwave Emission 



P^dv = ATTCabs{y)B^{T) 



(38) 



8.1. Candidate Materials 



behavior as in Fig. |2|; 



2 



100% of the Fe incorporated into small Fe304 grains; 



3 



5% of the Fe incorporated into bare Fe grains; 



4, 



5% of the Fe incorporated into small metallic inclusions. 
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For each case we assume a grain temperature T = 18K. IfV is the volume per H atom for a grain 
component, then its magnetic dipole contribution to the emissivity per H atom is 

ju ^ rigrV A'Khv'^ 1 9At2 .„„^ 

nH ~ UH c3 [eMhv/kT) - 1] {fxi + 2)2 + ^2 ' i i 

where we must remember that there may be additional electric dipole emission [see eq. (|l]-^)]. 

The magnetic dipole emissivity for each grain component is shown in Figure |^ also shown is 
the estimate for the "vibrational" or electric dipole emission, 

Sensitive studies of the microwave sky brightness have revealed microwave emission from 
interstellar matter; the emissivity per H nucleon has been deduced from the cross-correlation of the 
microwave sky brightness with far-infrared emission from dust grains, using measurements from 
the COBE DMR (Kogut et al. 1996), the ground-based Saskatoon experiment (de Oliveira-Costa 
et al. 1997), and OVRO (Leitch et al. 1997). These observational results are shown in Fig. |^ 

We see that a population of Fe grains or grains with Fe inclusions would be predicted to 
produce very strong emission near 70 GHz, with the peak occurring near the Frohlich resonance 
for Fe spheres. Even if the permeability of Fe has been estimated correctly, this peak would be 
broadened if, as expected, the Fe grains or inclusions are not all spherical, since the Frohlich 
resonance frequency depends upon the particle shape. The 90 GHz emission reported by Kogut 
et al. (1996) appears to limit the amount of interstellar Fe in metallic form to perhaps ^ 5% of 
the total Fe. In particular, this limits the fraction of interstellar Fe which can be present in Fe-Ni 
inclusions such as found in "GEMS" (Bradley 1994; Martin 1995). 



8.2. Hypothetical Materials 

We have estimated the emission expected from pure iron and magnetite, but the true state of 
Fe in interstellar grains is not known. It is possible that the Fe is concentrated in a form which 
is magnetic - not so strongly as pure iron, but more strongly than ferrimagnetic magnetite. This 
could, for example, be an Fe/Ni alloy with an appreciable concentration of O, H, Si, or other 
impurities. We now ask whether it is possible for magnetic dipole emission to account for the 
observed emission in the 14-90 GHz range. We approach this question by seeking to "customize" 
the magnetic properties, within the range of reasonable parameters, to see whether there are 
possible values which would lead to the observed emission. We make the following assumptions: 

1. The magnetic material is in (spontaneously-magnetized) single-domain grains, with a 
susceptibility x± perpendicular to the direction of spontaneous magnetization. 

2. The frequency dependence of the susceptibility is given by eq. (|26|). 

3. The low frequency susceptibility x±(0) should be somewhat smaller than the value 
X±(0) = 3.3 which we have estimated for pure iron. 
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10 102 
Frequency (GHz) 

Fig. 7. — Predicted thermal emission per H nucleon of interstellar dust for (1) paramagnetic silicate 
grains containing 100% of the Fe, Mg, and Si; (2) Fe304 grains containing 100% of the Fe; (3) 5% 
of the Fe in metallic iron, either as "bare" single-domain grains or as single-domain inclusions (see 
text) . Also shown is the estimated "vibrational" electric dipole emission from interstellar dust with 
^afe? ^ ^^'^ rotational emission from very small spinning grains as estimated for the "Cold 
Neutral Medium" by Draine &: Lazarian (1998a, b). Observed emissivities are inferred from data of 
Kogut et al. (1996; open diamonds), de Oliveira-Costa et al. (1997; filled diamonds), and Leitch et 
al. (1997; filled squares). 
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4. The spontaneous magnetization should be appreciably smaller than the value AirMg = 22kG 
for pure iron. As a result, the characteristic gyrofrequency loq/Itt = (e/mec)(2Ms/3) should 
be appreciably smaller than the ~ 20 GHz value estimated for Fe. 

5. We assume that nearly 100% of the interstellar Fe is in the hypothetical material X, and 
assume the Fe within the material X to contribute a mass density of 4 g cm^^ [thus the 
volume per H atom of the hypothetical material is Vx = 7.5 x 10^^^ cm'^/H - see eq. (|20|)]. 

Calculations are shown for 4 hypothetical materials, denoted XI, X2, X3, and X4, with x±(0) 
and ujQ values as shown in Fig. |8|. In Fig. ^ we see that in order to contribute an appreciable 
fraction of the observed 14-90 GHz emission, the hypothetical magnetic material must have 
X±(0) ~ 2 and ujq/2'k ~ 6 GHz. Note that material XI does not have a Frohlich resonance, and 
materials X2,X3,X4, with x±(0) = 2, have only a weak Frohlich resonance (see footnote ^). 

We conclude that if a large fraction of the interstellar Fe were incorporated into a moderately 
strong magnetic material with properties approximating those of our hypothetical material X4, 
then the thermal emission from these grains would approximately reproduce the observed thermal 
emission in the 14-90 GHz range. We stress that we are not arguing that such material exists 
- indeed, we have previously shown (Draine & Lazarian 1998a,b) that the observed emission 
seems likely to be due to rotational emission from a population of ultrasmall grains believed to be 
present for independent reasons. The point of the present discussion is to show that one cannot 
exclude the possibility that magnetic grains make an appreciable contribution to the emission in 
this frequency range. This issue should be resolved through observations of dark clouds (see §|To|). 

9. Polarization 

From Fig. |^ it is seen that if a small fraction of interstellar Fe is in iron grains, these could 
account for a substantial fraction of the diffuse emission in the 50 — 100 GHz region. Alternatively, 
we see in Fig. |^ that if a large fraction of the Fe is in a magnetic material with the properties of 
our hypothetical material "X4" , then a substantial fraction of the 14-90 GHz emission could be 
magnetic dipole emission from such grains. 

In recent years there has been significant progress toward understanding the alignment of 
interstellar grains. Analysis of the "crossover" phenomenon for grains subject to superthermal 
rotation now indicates that even ordinary paramagnetic relaxation may suffice to align larger 
(a ^ 10~^ cm) grains (Lazarian & Draine 1997), but radiative torques due to starlight appear 
to dominate those due to paramagnetic relaxation (Draine & Weingartner 1997). If smaller 
(a <^ 10^^ cm) strongly magnetic grains are present (e.g., iron or our hypothetical material X4), 
it appears likely that ferromagnetic relaxation would effectively align their angular momenta with 
the galactic magnetic field Bq, so that their "long" axes would tend to be perpendicular to Bq. 
If such grains are present, we then expect frequency-dependent polarization of the microwave 
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Frequency (GHz) 



Fig. 8. — Predicted thermal emission per H nucleon of interstellar dust for hypothetical Fe-rich 
materials X1-X4. In each case we assume that a fraction /x = 1 of interstellar Fe is in a substance 
with pxZx = 4gcm~^ [see eq. (pO])], with an assumed static susceptibility xo and characteristic 
gyrofrequency ujq = lO^wogs"^. Observed emission is as described in Fig. |^. Thermal emission 
from hypothetical material X4 could account for a substantial fraction of the observed 14-90 GHz 
emission. Also shown is the estimated "vibrational" electric dipole emission from interstellar dust 
with ) oc 
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emission from interstellar dust. 

To estimate the likely polarization, we assume the grains to be ellipsoids with semiaxes 
ffli < ^2 < ^3, spinning with the short axis di parallel to the angular momentum J (as expected 
for suprathermally rotating grains - see Purcell 1979). We assume "perfect" alignment of J with 
the galactic magnetic field Bq, with Bq perpendicular to the line-of-sight; the polarization of 
emitted radiation will be taken to be positive when the electric vector is perpendicular to Bq. For 
imperfect alignment, the polarization for perfect alignment should be multiplied by the "Rayleigh 
reduction factor" (Lee & Draine 1985) 

i?= (3/2)((cos^0) - 1/3) , (40) 

where is the angle between Bq and J. 

The predicted polarization depends on assumptions concerning the grain structure. We will 
consider two limiting cases. 



9.1. Single Domain Grains 

Suppose the grains which dominate the 10 — 100 GHz emision each consist of a single magnetic 
domain.^ We will assume the grains to have spontaneously magnetized with along the long 



axis 03 (this minimizes the magnetic energy). As discussed above in §7J., the permeability /i 
is anisotropic: jU = 1 for H || M^, and ^ = for H _L M^. Thus the degree of polarization 
P={Ie- Ih)/{Ie + h) where 



2 



£2 , €2 



|1 + L2(6-1)|2 |l + L3(e-l)|2 



+ (41) 

^ |l + Li(/xx-l)P ^ ' 



J €2 1 ^±,2 

h |l + Li(6-l)|2^2|l + L2(/i±-l)|2 ' ^ ^ 

The dielectric function of metallic iron can be approximated at low frequencies by a Drude model 
with ojQ = 0: 

- (13) 

where r ~ 3.8 x 10~^^ s and ujpT ~ 200. This approximately reproduces the tabulated values of 
e (Palik 1991). With this dielectric constant, the microwave emission is dominated by magnetic 
dipole radiation (|e| > 10^ for v < 100 GHz), so that < C^jf^ and the precise value of e is 

not critical. 

The resulting frequency-dependent polarization P is shown in Fig. ^ for different grain shapes. 
For both Fe or X4 grains, at high frequencies (z^ ^ 100 GHz for Fe, v ^ 25 GHz for X4) the 



^ Recall that Fe grains smaller than Oc ~ 3 x 10 ^ cm always consist of a single domain (Morrish 1980) . 
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polarization is large {P = 1/3) and positive (i.e., E perpendicular to the "short axis" of the grain). 
At these high frequencies \fi± — Ip <C 1 so that the values of the "shape factors" Li > L2 > L3 
are unimportant - the calculated polarization arises because the thermal fluctuations in the 
magnetization are perpendicular to 03; since 03 _L J the magnetic dipole emission tends to have 
H II J, and hence E _L J, for "positive" polarization. 

At lower frequencies 50 GHz for Fe, ^ 14 GHz for X4) |^ — 1| > 1 and shape effects now 
matter. There are still no magnetic fluctuations parallel to 03, but if oi is sufficiently small 
compared to 02, Li will be considerably larger than L2, and the magnetic dipole emission will 
tend to have H ± ai, hence E || ai, for "negative" polarization. This explains the behavior seen in 
Fig.i 

9.2. Grains With Magnetic Inclusions 

Suppose the grains to be ellipsoids containing isotropically-oriented magnetic inclusions with 
filling factor (j) ~ 0.03, as found for GEMS (see j |6.3| ). Since the grain material is now isotropic, the 
polarization is given by eq. (|8[jlO|). 

If the inclusions are metallic Fe, and contain N > 10^ atoms, we can use the permeability 
fj-{uj) shown in Fig. |5|. The resulting polarization for perfectly-aligned spinning grains is shown in 
Fig. |9|. Thermal emission from these grains is predominantly magnetic dipole radiation; since the 
grain material is isotropic]] 

the emitted radiation tends to have H perpendicular to the short axis di. However, since 
l/i — 1| <^ 1 (see Fig. |5|), the Lj-dependent term in eq. (^-|lO|) is only a small correction; hence 
I^K, Ij^K, (e2 + ^2) and the polarization P = {1^ — Ih)/{Ie + Ih) is small. 

10. Discussion 

Fe is the 9th most abundant element by number (after H, He, O, C, N, Ne, Mg, and Si), and 
5th by mass (after H, He, O, C), and nearly all of the interstellar Fe is in dust grains. It is therefore 
inevitable that some fraction of the interstellar grain material must be quite Fe-rich. We don't 
know what mineral form the Fe mainly resides in (a number of possibilities are given Table I), but 
the material will be at least paramagnetic, and quite possibly ferrimagnetic or ferromagnetic. 

It is therefore important to consider the emission of electromagnetic radiation as a result of 



^If the number of magnetic inclusions in a grain is small, the material within each grain will not, strictly speaking, 
be isotropic. However, if the orientations of individual inclusions are statistically independent of the overall grain 
shape, it is appropriate to approximate the grain material as isotropic in order to estimate the emission from the 
grain ensemble. 
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Frequency (GHz) 

Fig. 9. — Polarization of thermal emission from single-domain grains (upper panel) consisting of 
either metallic Fe (solid lines) or the hypothetical material X4 (broken lines), assuming perfect 
(spinning) alignment (see text), Single-domain prolate spheroids (oi : a2 : 03 :: 1 : 1 :> 1) have 
P = 1/3 for u ^ 100 GHz. The lower panel shows the polarization of the thermal emission from 
perfectly aligned ellipsoids composed of material with metallic Fe inclusions with volume filling 
factor (f) = 0.02. Results are shown for various values of the axial ratios ai :o2 :a3. 
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thermal fluctuations in the magnetization of the grain material. These fluctuations take place at 
microwave frequencies and below (the magnetization is constant on timescales short compared 
to the time for electrons to precess in the internal magnetic fields within the material; these 
precession frequencies are at most ~ 20 GHz (the value for metallic iron). As a result, thermal 
magnetic dipole emission from interstellar grains ~ 100 GHz may be be stronger than the thermal 
electric dipole "vibrational" emission (arising from thermal fluctuation is in the charge distribution 
within the grain) at frequencies <^ 100 GHz.^ We show how this magnetic dipole emission can be 
calculated, but the emission requires knowledge of the magnetic permeability at these frequencies. 
We have estimated this permeability for various materials of interest, including metallic iron and 
magnetite (Fe304), but have been required to extrapolate to the high frequencies of interest. 
It would be of great value to have direct measurements of the magnetic properties of small Fe 
particles at 50-100 GHz. 

According to our estimates, metallic iron would produce such strong emission near ~ 70 GHz 
that not more than ~ 5% of the Fe can be in the form of pure metallic iron in order not to exceed 
the observed emission at these frequencies. This limits the fraction of interstellar grains with Fe-Ni 
inclusions such as in the GEMS found by Bradley (1994). 

Magnetite, on the other hand, would radiate only a fraction of the power emitted by 
interstellar grains in the 14-90 GHz region. However, it is possible that interstellar Fe could 
be present in some Fe-rich substance (e.g., an Fe/Ni alloy with an appreciable concentration of 
Mg,Si,0,H impurities) with magnetic properties intermediate between those of metallic iron and 
magnetite, and we therefore cannot rule out the possibility that the bulk of the observed emission 
in the 14-90 GHz range could be thermal magnetic dipole radiation. It would be of great value to 
have laboratory measurements on various plausible candidate materials (e.g., amorphous olivine, 
or Fe-containing alloys) at these frequencies. 

Rotational emission from spinning dust grains has previously been proposed as the explanation 
for the observed 15-90 GHz emission from dust in diffuse clouds (DL98a,b). There are two ways 
in which this new emission mechanism - magnetic dipole radiation from thermal fluctuations in 
grain magnetization - can be distinguished from the hypothesized rotational emission: 

• The rotational emission requires ultrasmall grains. Since we have reason to believe that 
ultrasmall grains are depleted in dense regions, we would then expect the rotational emission 
to be relatively weak in dense gas. On the contrary, thermal magnetic dipole emission from 
dust should be largely unaffected by coagulation of the dust grains. Therefore, observations 
of (or upper limits on) the 10 — 60 GHz emission from dense clouds can be distinguish 
between rotational emission from ultrasmall grains or thermal magnetic dipole emisson from 
magnetic grain materials. 



Rotational electric dipole radiation from rapidly-rotating ultrasmall grains is also important at these frequencies 
(Draine & Lazarian 1998a, 1998b). 
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• If produced largely by single-domain grains, the magnetic dipole emission potentially 
has a complex and strong polarization signature, which could be very different from the 
polarization expected from spinning dust grains which have been partially-aligned by 
magnetic dissipation (Lazarian & Draine 1998). 

11. Summary 

The principal results of this paper are as follows: 

1. Formulae are presented for electric and magnetic dipole absorption cross sections for 
homogeneous spheres or ellipsoids of material with dielectric function e and magnetic 
permeability fi. 

2. The Kramers-Kronig relation relating the total grain volume to /q°° CextdX is generalized 
to include the case of magnetic grains. The resulting eq. (|l^) shows that if an appreciable 
fraction of the grain material is magnetic, then magnetic dipole effects must dominate the 
extinction cross section at the frequencies ly ^ 30 GHz where plausible grain materials have 
a magnetic response. 

3. Because most of the Fe is in solid form in the interstellar medium, it is expected that 
some fraction of the interstellar grain population must be appreciably magnetic - either 
paramagnetic, superparamagnetic, ferrimagnetic, or ferromagnetic. The frequency-dependent 
magnetic properties of candidate materials are estimated. 

4. The magnetic dipole emission from either paramagnetic grains or magnetite (Fe304) grains 
is well below the expected electric dipole emission due to vibrational modes at 100 GHz, 
or rotational emission from very small grains at ^ 100 GHz. 

5. Fe grains or inclusions will have the magnetic analogue of a Frohlich resonance at ~ 70 GHz, 
where the magnetic dipole absorption cross section peaks because //i = —2. This will result 
in strong magnetic dipole absorption from ferromagnetic iron particles near ~ 70 GHz. 

6. If our estimate for the magnetic susceptibility of single-domain Fe is correct, not more than 
~ 5% of interstellar Fe can be in the form of metallic iron; otherwise the thermal magnetic 
dipole emission from interstellar dust at 90 GHz would greatly exceed the emission from dust 
measured by Kogut et al. (1996). 

7. If a substantial fraction of interstellar Fe is in a moderately strong magnetic material with 
the properties resembling our hypothetical material "X4" in Fig. ^, then thermal emission 
from such grains could account for a substantial fraction of the observed 14-90 GHz emission. 
At this time this possibility cannot be excluded, even though the 14-90 GHz emission seems 
most likely due to spinning ultrasmall dust grains (Draine & Lazarian 1998a, b). 
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8. If nonspherical single-domain ferromagnetic grains are present, the magnetic dipole emission 
will be polarized. Such grains are expected to be spontaneously magnetized along the 
"long" axis. The polarization is expected to depend strongly on frequency, with "normal" 
polarization (electric vector perpendicular to Bq) at frequencies ^ 100 GHz, but with a 
decrease, and perhaps even reversal (depending upon grain shape), of the polarization for 
^ 30 GHz. 
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A. Microwave Absorption due to Stark-EfFect Splitting 

The electric fields within a solid - whether crystalline or amorphous - can split the magnetic 
sub-levels of paramagnetic ions in the solid. Kittel & Luttinger (1948) noted that while these 
splittings usually correspond to frequencies ~ 10^^ Hz, under some circumstances the splittings 
correspond to microwave frequencies. To estimate the magnitude of the associated absorption, 
suppose that the solid contain a density dn of species with energy levels split by frequencies in the 
interval duj, and let \fiui\'^ be the magnetic dipole matrix element between the upper and lower 
states, giving a spontaneous decay rate 

Aui = ^If^uil'' ■ (Al) 
If u; ^ kT/h, then d{ni — Uy) = {hui/2kT)dn. The absorption coefficient is 

27r^ fhu;\ dn l^nZp 

" = ^lfcrjdi^^ ' 

corresponding to an imaginary component of the magnetic susceptibility 

c 27r2 / dn 
A*2 = -a = -hrFT Ti (A3) 




2.7x10- ^ . (A4) 



where //s is the Bohr magneton. The total atomic density ~ 10^^ cm~^ in the solid, so 
dn/dlnu = 10^^ cm~^ would correspond to 10% of the atoms having resonances within a frequency 
interval A In a; = 1. Such a high value would probably only occur if the Stark-effect splittings are 
concentrated near a few frequencies. 
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B. Susceptibility of Single-Domain Ferromagnetic Particles 

Consider a spherical particle consisting of a single domain. At temperatures T ^ Tc/2, where 
Tc is the Curie temperature, we may take the spontaneous magnetization to be approximately 
equal to the saturation magnetization Ms- Suppose the crystal has cubic symmetry. The free 
energy is a function of the direction of the magnetization. If ai, 02, and are the direction 
cosines of the magnetization relative to the crystal axes, the "anisotropy" free energy is 

Fk = Ki{a\a\ + a\a\ + a\oi\) + K2a\a\a\ (Bl) 

where Ki and K2 are the anisotropy constants (Morrish 1980). 



B.l. i^i > 0, > 

If Ki > and K2 > 0, then the "easy" directions are along the cubic axes (e.g., (100)). If 
we now apply a weak magnetic field H± perpendicular to the magnetization, the magnetization 
will deflect by an angle 9. For weak fields, ^ ^ 1 and the anisotropy energy and magnetic energy 
become 

Fk = KiO^ + 0{e^) (B2) 

Fh = -M • H = M,H±e + 0(6|2) . (B3) 

Minimizing [Fk + Fh)-, we find 9 = MsH±/2Ki. Thus, the transverse susceptibility x± = Ms6/H± 
is 

For Fe, Ms = 1750 G, and Ki = 4.6 x lO^ergscm'^ (Morrish 1980), thus x±(0) = 3.3 and 
/^±(0) = 1 + 47rx±(0) = 43. This is significantly smaller than the value ;u(0) ~ 150 characterizing 
bulk (multidomain) Fe. 



B.2. Ki <0, K2 <0 

When Ki < 0, K2 < 0, the "easy" directions are along the diagonals (e.g., (111)). The 
anisotropy energy is now 

with Fh again stih by eq.(p3|). Again minimizing {Fk + Fh) we find 9 = -9MsH±/[12Ki + AK2], 
and 

^^(0) « ^ (B6) 

' I2K1+AK2 ^ ' 
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Maghemite (7Fe203) is cubic with 47rAfs = 4780 G (Dunlop & Ozdemir 1997). The magnetic 
properties do not appear to have been measured at low temperatures; at room temperature, 
Ki = —2.5 X lO^ergscm"^ (Valstyn, Hanton, & Morrish 1962), from which we estimate 
X±(0) ~ 0.43. As discussed in §^ above, a shghtly larger value of x±(0) ~ 0.6 is in agreement 
with measurements by Valstyn et al. (1962) at GHz frequencies. 

B.3. Oscillating Fields 

The above discussion applies to static fields. For oscillating fields, we estimate x±('^) using 
eq. (p6D, with ujq = (e/mec)(47rMs/3), the precession frequency of an electron in the internal field 
47rM,/3. 
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